from qutip import *
import numpy as np
import matplotlib.pyplot as plt
sys.path.append('D:\Dropbox\Dokumente\PI3\Code\Python\quantum-tools')
import qutip_enhanced
bloch = Bloch3d()
qte = qutip_enhanced.QutipEnhanced()
import nv_hamilton as nvh

nvh = nvh.NVHam(magnet_field = {'z': 0.55}, n_type = None, electron_levels = [1,2])
# nvh.add_spin(nvh.hft_13c_dd(location={'rho': 1e-9, 'elev': np.pi/3., 'azim': 0}), nvh.h_13c(), [0, 1])
print nvh.h_nv
# print (a[0] - a[1]) - (a[2] - a[3])